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Abstract 



We consider finite volume (or equivalently, finite temperature) expectation values of 
local operators in integrable quantum field theories using a combination of numerical 
and analytical approaches. It is shown that the truncated conformal space approach, 
when supplemented with a recently proposed renormalization group, can be sufficiently 
extended to the low-energy regime that it can be matched with high precision by the 
low-temperature expansion proposed by Leclair and Mussardo. Besides verifying the 
consistency of the two descriptions, their combination leads to an evaluation of expectation 
values which is valid to a very high precision for all volume/temperature scales. As a side 
result of the investigation, we also discuss some unexpected singularities in the framework 
recently proposed by Pozsgay and Takacs for the description of matrix elements of local 
operators in finite volume, and show that while some of these singularities are resolved 
by the inclusion of the class of exponential finite size corrections known as /i-terms, these 
latter corrections themselves lead to the appearance of new singularities. We point out 
that a fully consistent description of finite volume matrix elements is expected to be free of 
singularities, and therefore a more complete and systematic understanding of exponential 
finite size corrections is necessary. 



1 Introduction 

Finite temperature expectation values of local observables are important for numerous appli- 
cations of quantum field theory. In 1 + 1 dimensional integrable quantum field theories, they 
have been the subject of intensive studies recently. In [I], Leclair and Mussardo conjectured 
a low-temperature expansion which uses two ingredients: the thermodynamic Bethe Ansatz 
[2j as applied to integrable quantum field theories [3], and the form factor bootstrap program 
[H E| • Although a partial proof of the series was given by Saleur in [7] , it was not until the 
recent introduction of the finite volume form factor framework by Pozsgay and Takacs [U [9] 
that a systematic construction of the series itself was made possible. The finite volume form 
factor formalism even led to an explicit derivation and further generalization of the series to 
all orders [10]. The series itself has found applications in the investigations of one-dimensional 
quantum gases |11[ [T2| 113] , The finite volume form factor formalism, on the other hand, was 
extended to two-point functions at finite temperature |14( [To"] and a method to determine the 
expansion to any given order was developed in |16[ IT7] . 

The aim of the present work is to address some open issues, regarding both the finite volume 
form factor formalism and the determination of thermal one-point functions. Regarding the 
latter, we make use of the equivalence between a system at finite temperature T and in a finite 
volume L = 1/T, we consider the one-point function in finite volume, as determined using the 
following methods: 

• the Leclair-Mussardo series; 

• the thermodynamic Bethe Ansatz (TBA), which is equivalent to the Leclair-Mussardo 
(LM) series applied to the trace of the energy-momentum tensor; 

• and the truncated conformal space approach (TCSA), introduced in [18j . 

The LM series is a low-energy expansion, valid for large values of L, while the TCSA is a non- 
perturbative extension of the ultraviolet perturbation theory, which is valid for small values 
of L. Therefore care must be taken to find a regime where these approaches can be matched. 
This is facilitated by applying a recently introduced TCSA renormalization group method 
|19[ |20| I2T] to the evaluation of one-point functions, in order to improve the precision of their 
TCSA evaluation for larger values of the volume. 

In addition to the thermal/finite volume expectations values, this work also discusses a 
separate, but related issue that we found when testing for the consistency of the TCSA and 
the bootstrap form factors. Namely, due to singular behaviour in the density of states, the 
finite volume form factors predicted by the formulae in [8, 9 J show some singularities. It turns 
out that these can be resolved by including some of the exponential finite size corrections, the 
so-called /x-terms along the lines of |22| . 

For concrete calculations we chose the T2 model which is an integrable perturbation of the 
minimal conformal field theory A4.2.1 by its relevant perturbation $1,3. This was motivated 
by the fact that this is one of the simplest theories in which there is an additional non-trivial 
primary field besides the perturbing one, so this provides a non-trivial application of the LM 
series in the sense that the series for this operator is not a consequence of the TBA. 

The outline of the paper is as follows. In section [2] we briefly recall relevant facts about 
the LM series and TBA. In section [3] we give a derivation of the cut-off dependence of expec- 
tation values in TCSA, since we need the explicit exponents for our calculations. While the 



extrapolation has already been used on at least two occasions |20[ [23] , the theory behind the 
exponents has never been exposed and our derivation fills this gap. Section [J] is devoted to 
the comparison of the one-point function obtained from the TCSA on the one hand, and the 
LM series (or, equivalently, from the TBA in the case of the trace of the energy-momentum 
tensor) on the other. In section [5] we digress a little to discuss the issue of singularities in the 
finite volume form factor formalism and their resolution by the inclusion of //-term corrections. 
Section [6] is reserved for the conclusions. There are two appendices: Appendix [A] gives the 
relevant details about TCSA and its application to the T2 model, while Appendix [B] is a brief 
summary of its scattering theory and form factors. 

2 One-point functions at finite temperature 

2.1 The Leclair-Mussardo series 

2.1.1 The LM series itself 

The conjecture for the finite volume expectation value of a local operator by LeClair and 
Mussardo [lj takes the following form in the T 2 model: 



y, _L_ r° d0i d^d^i d ^TT 1 TT l 

{ )l ^ n!m! Loo 2vr " ' 2vr 2vr " ' 2vr ±\ 1 + e £ i(^) 1X 1 + p ^(§j) 

n,m=0 i=l ]=1 ± ~r c 
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where 9i, 9{ are rapidities of the particles with mass mi and 777-2. The e a are the pseudo-energy 
functions that satisfy the TBA integral equation 

— ip ab {6-e')\og(l + e- £ ^) (2.2) 

with ip a b the derivative of the phase shift of the S-matrix, and F^ 2m c is the connected form 
factor defined by 
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where the finite part is taken by expanding the e-dependence of the form factor expression for 
small values of the e-s and omitting all terms which blow up when any of the et is taken to 
(including terms that depend on the ratios of the e/%) pQ. 

2.1.2 Evidence for the LM series: TBA and finite volume form factors 

Leclair and Mussardo showed that for the trace of energy momentum tensor their series coin- 
cides with the result derived from the TBA by Zamolodchikov [3], which we briefly summarize 
below specialized to the T2 model. 



Since Ti is the perturbation of M.2,1 minimal model by $13 operator, the trace of the 
stress energy tensor is given b>tj 



9 = T£ = 27v\(2h 1>3 - 2)$ lj3 (2.4) 

where A is the coupling constant and his the conformal weight of the operator. The ground 
state energy in finite volume is given by 
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where the e a solve fj2.2|) . The vacuum expectation value of can be expressed by the following 
derivative of Eq (L) 
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where the derivatives of the TBA pseudo-energy solve the following linear equations 
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By expanding the solution to these equations for large L, one obtains the series (|2.ip [T]. 
Extending the argument, Saleur proved the equivalence of the TBA and the LM series for 
one-point functions of conserved charge densities [7] . 

For general operators where a TBA construction is not known, the LM series was demon- 
strated to be true up to (and including) three particle terms in [9j, and a proof to all orders 
was presented in |10| . We remark that these results were proven using the finite volume form 
factor formalism introduced in [HI E]? and the proofs relies on the expression of finite volume 
diagonal matrix elements (up to exponential corrections in the volume), which is still only a 
conjecture, albeit well-supported both by analytical arguments 0[9] and numerical data [5]. 

3 One-point functions from TCSA using RG extrapolation 

3.1 Cut-off dependence of the vacuum expectation value 

For a numerical evaluation of the expectation values, we used the truncated conformal space 
approach (TCSA) developed by Yurov and Zamolodchikov in |18| ; for the technical details we 
refer the reader to appendix lAl 



lr The stress energy tensor is normalized following the conventions of Zamolodchikov's paper [3]. 



In order to test the exponentially decaying finite volume corrections predicted by the LM 
series, we need to improve the accuracy of the TCSA. This can be achieved by extrapolating 
in the cut-off (truncation level) dependence using TCSA RG, which was first applied in |20j, 
and later also in |23j . However, the way to derive the necessary exponents has not yet been 
presented; therefore here we give the necessary details below so that we can apply it to our 
case, and also in order to facilitate future applications of the method. 

The derivation of the cut-off dependence below uses the general ideas presented in [21 J 
where a similar method was used to obtain the running coupling constants. We parametrize 
the cut-off e cu t by the left descendant quantum number of the highest vectors included in the 
truncated space. In the T2 model, this is the same for all modules in the Hilbert space as all 
the differences between the primary dimensions are smaller than 2. The vacuum expectation 
value of the operator O, when cutting TCSA at level n is: 

E{n) - Q{n) 



M{n) 



Q{n) = (min| P n O(0,0)P n e- X f d2zP " Viz ' s)Pn \mm) 

where z = r — ix, V is the perturbing operator, which has conformal spin zero (i.e. hy = hy) 
to ensure translational invariance; to simplify the exposition, we also assume that ho = ho- 
All operator products are assumed to be time ordered, P n is the projector on states with 
descendant numbers less than or equal to n and |min) is the lowest weight primary state, 
whose perturbation becomes the vacuum of the massive theory (for T2 this is the state created 
by $1,3) . The normalization factor is 

M{n) = (min| e -*fd 2 zP n V(z,z)P n \ min j 

Our aim is to calculate the dependence on n, i.e. the difference Q(n) — Q(n— 1). To simplify the 
subsequent formulae, let us introduce the following notation for the past and future evolution 
operators corresponding to the cut-off Hamiltonian: 



4 n) (A) = Texp (-X J d 2 z'P n V(z',z')P n ] 



so we can make time ordering explicit 

Q(n) = (min|[/| n) (A)P n O(0,0)P n C/ W (A)|min) 

We assume that n is large, and so the states at descendant level n have energy much higher 

than the mass gap m: 
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When this condition is satisfied, we can evaluate their contribution to first order in A. Defining 
the projector on the level n descendants as P n = P n — P n -\ and denoting 

A n V(z, z) = P n -iV{z, z)P n + P n V(z, z)P n - X + P n V(z, z)P n 



the r < time-evolution operator can be expanded as 
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and similarly for IfY 1 (A). This yields 
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where we used that as a consequence of 
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P n commutes with LQ? (A), and also the relations 



P n |min) = 
P n _i | min) = 

Exploiting invariance under time reflection t ^ —t 



f/i n_1) (A) \ |min) 



mm) 



(n-1), 



Q(n)-Q(n-1) = -2A / d 2 z (min| ?7J n ~ ij (A) O(0,0)P n V(z,z) U { _~ l) (X) |min) + 0(A 



r(»-i)/ 



T<0 



Applying the exponential map (jA.ip and writing w = re 
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where we used that to leading order one can take n — > oo in the time evolution exponentials. 
Using the operator product expansion (OPE) in conformal field theory 
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where A runs over the set of scaling fields and the C are the operator product structure 
constants. 

As described in Appendix lAl due to the fact that the perturbing operator satisfied hy = 
hy, the TCSA Hilbert space decomposes into sectors with a given value of the conformal 
spin Lq — Lq, which is in fact a consequence of translational invariance. Since the vacuum 
state is translationally invariant it has zero conformal spin. Therefore all intermediate states 
contributing to the matrix element have spin zero, and so the operators A contributing to the 
cutoff dependence also have zero conformal spin, i.e. Ha = ^4 as well. Substituting (|3.3p into 
(|3.2p and using the expansions 

1 ^ T(a + n) 



r(a)T( 



(l- w )« ^r(a)r(n + l) 
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the (n, n) term just gives the contribution to the OV two-point function from intermediate 
states at level (n,n). Due to the angular integral the only contributing terms are those in 
which n = n, which is another manifestation of translational invariance. The exponentials 
ensure that the contribution of a given operator A in the OPE will be proportional to the 
exact finite volume vacuum expectation value in the perturbed theory 

(min| U { °°\x)A(w,w)U { _ QO \\) |min) 

w=w=l 

This captures the non-perturbative infrared physics, while perturbation theory applies in the 
UV regime and describes the dependence on n (provided that n is large enough). 

Therefore a given operator A contributes to Q(n) — Q(n — 1) a term proportional to 

1 / r(q A + n) \ 2 2n-l+2h vd = 1 ( r ( a A + n) x ' 

, \T{a A )T(n + l)J ' " 2(n + h v ) \T(a A )T{n + 1) 

where 

«A = ho + hy — Ha 

and all factors independent of n were omitted. To evaluate the leading behaviour for large n, 
one can use Stirling's formula 



T(x + 1) « V2^x x e~ x (1 + 0(l/x)) 

and the relation 

(n + a) (n+a) w n n+a e a (1 + 0(1/ n)) 

which results in 

Q{n) - Q( n - 1) oc ^ K A n 2aA ~ 3 (1 + 0(l/n)) 

where iC A are n-independent coefficients. For large n, this gives the differential equation^! 
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2 We remark that the difference equation is also easy to solve and the difference from the differential equation 
is numerically important for small n. On the other hand, for the relatively high values of the cutoff the difference 
between the two is insignificant. 



Solving the equation with the initial condition at n = oo, the result is 

Q(n) = Q(oo) + Y, K A n 2aA ~ 2 (1 + 0(l/n)) + 0(A 2 ) 
A 
The leading corrections in the descendant level are of order 1/n (as indicated), arising from 
omitting hy besides n, and also from the approximations used in Stirling's formula. 

For the complete expectation value E(n) one must also compute the running of the nor- 
malization factor JV"(n). This is the same as Q(n) for the choice when O is the identity, and so 
A = V is the only operator in the OPE. Therefore the relevant exponent a is ho + hy — Ka = 
and there is no w, w dependence in the OPE, and so J\f(n) does not run to first order in A: 

^ = 0(X 2 ) 
an 

Supposing that the expectation value we compute is convergent in TCSA (i.e. 2a A — 2 < 

for all A), the end result for the matrix element is: 

E(n) = E{oo) + J2 K A n 2aA ' 2 (1 + 0(l/n)) + 0(A 2 ) (3.4) 

A 

The issue now is whether this calculation is self-consistent, i.e. whether the 0(A 2 ) terms are 
smaller than the O(A) terms we evaluated. Since they are rather complicated to evaluate, 
we rather use a "proxy" condition which is to require that the first order correction itself is 
small when compared to E(oo). One can neglect the exponentially small difference between 
the infinite and finite volume expectation values of O: 

25(oo) - (0) x ~ O (e~ mL ) 

and assuming that A is a primary field (so it transforms homogeneously under the exponential 
mapping) we get 

(T \ 2/lA 
— (A) x + O (e- mL ) 

2ir J 

in terms of the infinity volume vacuum expectation value (A) x . Up to some dimensionless 
numerical coefficients, this leads to 

Simple dimensional analysis gives 

Aocm 2 " 2 ^ , (A) x ocm 2 - 2AA , (0) x oc m 2 ~ 2A o 
up to numerical coefficients. This leads to 
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which is the same condition as ([3.ip . This is indeed consistent with what we have observed 
in all of our practical calculations performed so far: the first order extrapolation breaks down 
in large enough volume (the precise value of which depends on the model considered); and 
the higher the value of the cut-off, the larger is the critical value of the volume where the 
break-down occurs. 



3.2 Remarks 

First of all note that to leading order all the calculation is unchanged if we consider a general 
matrix element instead of a vacuum expectation value, as long as the cut-off is chosen so that 
the states between which the matrix element is taken are included in the truncated Hilbert 
space. The only difference is that instead of the vacuum matrix element of operator A one 
obtains the appropriate excited state finite- volume matrix element as coefficient. Therefore, 
while the exponent of the cutoff dependence is universal, its coefficient depends on the matrix 
element considered. 

If the operator is such that its matrix elements are not convergent with increasing cutoff, 
the operator must be renormalized by subtracting the divergent parts. These are always 
proportional to operators with lower scaling dimensions, as long as the perturbation is relevant 
(ie. the theory is renormalizable) . This is true because the condition for divergence is 

2a A - 2 > 

which leads to 

2(ho + h v - h A ) - 2 > 

and since in a renormalizable theory the perturbing operator is relevant (hy < 1), a divergent 
term can only be obtained if for some A 

h A < h a 

In this case, a renormalized operator [O] must be defined by subtracting from O a counter term 
proportional to A. This is consistent with the general theorems of renormalization derived in 
the framework of the standard Feynman perturbation expansion in renormalizable quantum 
field theory. 

Finally, note that the cut-off dependence can also be obtained by a simple scaling argument, 
by examining the operator product expansion: 
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A 
The energy dependence of the coefficient is given by simple dimensional analysis as 

C A oc [energy] " 2+2aA , a A = h + h v -h A 

The typical energy of (zero-spin) states at descendant level n is 

Aim 
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If n is high, then this energy is much larger than the mass scale m, so by usual rules of 
ultraviolet perturbation theory the dependence on the mass scale can be neglected to leading 
order, which means that the above energy is the only scale in the calculation. Therefore the 
contribution from level n states to leading order takes the form 



C A oc 



Aim 



-2+2q 



A 



L 
which also gives the correct volume dependence as well. 



In TCSA everything is measured in units of the mass scale m, so the relevant dimensionless 
number which characterizes the magnitude of the first order correction is 

n 2-2a A f^n\- 2+2aA 
c A = C A m 2 2aA x I — J 

in agreement with (|3.5p . 



3.3 Some particular models 

3.3.1 Ising model in a magnetic field 

This is the model with the famous E% scattering theory obtained by Zamolodchikov |24j . The 
Hamiltonian is given by 

H = H* + h I dxa 



where a is the spin-field with h = h = 1/16. For the matrix elements of a itself, the most 
relevant operators occurring in aa are 

A = \ : h A = 
A = e : h A = l/2 

and so the corresponding leading cut-off dependencies are 

n" 7 / 4 , n~ n / 4 

where the second exponential is of the same order as the (neglected) 1/n corrections to the 
first. In fact, an extrapolation using the first exponent gives an excellent result for the vacuum- 
two-particle matrix element, as demonstrated in [20J. 

3.3.2 Sine-Gordon model 

This is a perturbation of a massless free boson (c = 1 conformal field theory) with the Hamil- 
tonian 

H = I dx]- Udt$f + (d x $) 2 ) -A f dx : cos /3$ : 

and so in this case V =: cos/33? :. 

At first sight our derivation of the cut-off dependence is not really applicable to the sine- 
Gordon model. As the spectrum of primary fields is not bounded from above, with increasing 
cut-off e cu t new Fock modules enter the Hilbert space |25j . Therefore there is no simple relation 
between the cut-off and the level of the highest descendants in the Hilbert space. Nevertheless, 
it turns out that in the case discussed below only two operators contribute in the operator 
product, and therefore there is a one-to-one relation between e cu t and the highest descendant 
level n included in the two contributing Fock modules, which is enough for our derivation to 
apply. 



Let us investigate the operator 
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We have 

ho = hv = h = — 

Using the operator product of exponential fields 

. e ia$(z,z) .. e ip§(w,w) ._ \ z _ w \-t£ ■ e ia<5>(z,z)+il3<5>(w,w) . 



j leading operators in the OV OPE are: 




A = 1 : h A = 


= 


A = d$8$ : h A = 


= 1 


A = e*P* : h A = 


= 4h 


ich lead to the following cut-off dependence: 




n -2+4h n -4+4h 


, n~ 2 ~ Ah 



Note that for 8h > 1, the 1/n corrections to the first term have a larger exponent than the 
third. 

The extrapolation with only the leading term, when combined with a numerical renormal- 
ization group, gives excellent results in the attractive regime, as demonstrated in [23J. 

3.3.3 Scaling Lee- Yang model 

This is the model for which TCSA was originally developed |18| . In this case, there is a single 
non-trivial primary field <J>. Choosing <J> as both the perturbing and measured operator: 

o = v = $ 

one has 

ho = hy = h = — - 
5 

The operators occurring in the OV are 

A = $ : h A = -- => o-a = h 
5 

A = l : h A = 0=> a A = 2h 
and so the cut-off dependence is given by the terms 

When applying this cut-off dependence to extrapolate TCSA measurements of the finite vol- 
ume vacuum expectation value to infinite cut-off, it turned out that TCSA data are so detailed 
and accurate that we could also fit the subleading exponents 

n -17/5 ; n -19/5 

resulting from the 1/n corrections. These fits to the cutoff dependence are very accurate, 
and the resulting extrapolating matrix elements are shown in figure 13.11 These results also 
demonstrates that the RG improvement is really necessary to guarantee a proper matching 
between TCSA and the Leclair-Mussardo series (which is equivalent to the TBA in this partic- 
ular case). Note that, as predicted by the argument in Section T3. 21 the extrapolation becomes 
inaccurate in large volume. One expects that in this case it is necessary to evaluate the cutoff 
dependence to second order in perturbation theory, which woud require integrating a four 
point function in the conformal field theory, and is out of scope of the present work. 

11 
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Figure 3.1: Extrapolation of the dimensionless vacuum expectation value im 2 ' 5 ($1,3) in the 
scaling Lee- Yang model, compared to the exact value computed from TBA. 

3.3.4 T 2 model 

The model itself is described in Appendices [A] and |Bj Let us first take the measured operator 
to be 

O = $i i2 
Since V = $1,3, the operators occurring in the OV are 



A = $1,2 : h A = --=>a A 
A = $i >4 : h A = -■= => a A 
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and so the cut-off dependence is given by 
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For the case 



O = $1,3 
the operators occurring in the OV are 

A = 1 : h A = => a A 
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Figure 3.2: Extrapolation of the dimensionless vacuum expectation value m 6 ' 7 ($1,3) in the 
T2 model, compared to the exact value computed from TBA. 



and so the cut-off dependence is 



n 



-20/7 



?? 



-26/7 



The result of the extrapolation is for the perturbing operator $13 is shown in figure [3~2| again, 
the RG improvement is really necessary to guarantee a proper agreement. 



4 Numerical comparison 

Here we aim to compare the Leclair-Mussardo series to the RG-extrapolated TCSA results, 
for the operators $13 and $1,2- The expectation value of the $13 operator in the T2 model is 
related to the trace of the stress-energy tensor (|2.4|l , so its expectation value can be evaluated 
from TBA and it does not give us an independent validation for the LM series. However, the 
error in the numerical solution of the TBA equation (|4.1|) can be made very small (in our 
case it was of the order 10 -14 ). We used this to check our procedures for evaluation of the 
LM series, and found excellent agreement. Another application is to estimate the accuracy 
of the RG-extrapolated TCSA. Table 14.11 shows the difference of the expectation value for 
the $i ; 3 operator between the TBA prediction and RG-extrapolated TCSA evaluation in 
units of m\ for several values of the dimensionless volume parameter, I = m\L. The results 
show the efficiency of the RG-extrapolated TCSA. We also showed the difference between the 
RG-extrapolated TCSA value and the one evaluated at the highest value of the cutoff; as 
exemplified by the data ate I = 15, the extrapolation results of an improvement of accuracy 
of almost two orders of magnitude for volumes I > 10. 
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Table 4.1: The difference Sfis = m-/ ((^i,3)tba — ^ 1 > 3 )tcsa) between the RG-extrapolated 
TCSA evaluation and the TBA prediction for the VEV of $1,3. 



Now we turn to the evaluation of the LM-series for the $1^ operator. At infinite volume 
only the infinite volume VEV ()B.4p contributes to the LM-series, while for small volumes 
almost all terms are relevant. Their contributions are proportional to ~ e~^ , where n is the 
sum of the masses of the particles entering a given term. For every contribution there's a 
characteristic volume where it becomes relevant. We calculated the first seven finite volume 
correction terms: 1, 2, 11, 12, 111, 22 and 112, where the number of 1 and 2 means the 
number of particle with mass m\ and m.2 participating in the given contributions. Table 
4.21 shows the value of the RG-extrapolated TCSA evaluation and the LM-series with the 
mentioned contributions for $1,2, while Table 1131 shows the difference of the LM-series from 
the RG-extrapolated TCSA data while adding more and more contributions. 

For I > 1, there is a steady improvement as more and more terms are added, and it is 
clear from the table that contributions from higher states are switched on at progressively 
lower volumes. Comparing the columns labeled +111 and +12, it can be seen that adding the 
111 contribution makes the agreement of the LM series with TCSA worse; the reason is that 
the term 22 is of the same order of magnitude as the contribution from 111, so consistency 
requires adding them to the series together. Indeed, the deviation reported in column +22 is 
smaller than the one in column+12. 
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Table 4.2: The values for im 1 
Leclair-Mussardo series. 
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(^1,2) from the RG-extrapolated TCSA evaluation and the 



For I < 1, the series does not converge very well; indeed for I = 0.5, there is no sign 
of any convergence. However, for such small values of the volume higher terms of the LM 
series would still be significant. Indeed, the leading corrections following the 112 term are the 
contributions! Ill and 122, which can be estimated to be of order 



-4miL 



and 



-(mi+2m,2)L 



which at I = 0.5 give approximately 14% and 12%, respectively. This agrees well with the 
magnitude of the deviation at / = 0.5, as can also be seen from table l4~2l 

For volumes I < 3.5 the estimated TCSA error is smaller than 10 -6 , and so the deviation 
between TCSA and the LM series is dominated by the higher corrections to the LM series. 
However, the calculation of higher contributions becomes progressively slower as the number 
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((^i,2)lm — (^i> 2 )tcsa) between the RG-extrapolated 
TCSA evaluation and the Leclair-Mussardo series, depending on the multi-particle contribu- 
tions included in the latter. 



of particles to include increases (122 can be evaluated using a 10-particle form factor, cf. the 
remark at the end of Appendix IB]) . Similarly to the case of 111 and 22, contributions 1111 and 
122 are roughly of the same order and so they must be added to the series together. Treating 
the 10-particle form factor numerically proved to be rather difficult, so we have not evaluated 
the contribution 122. 

However, there is a way to verify the consistency of the above considerations further. We 
fitted the deviation in the last column of Table EOl the difference from the TCSA by the Ansatz 



ae 



-4miL 



+ be 



-(mi+2m2)L 



dictated by the particle content of the states 1111 and 122. As shown in figure I4~TI the fit 
was very successful. In addition, the parameters a and b turned out to be of the same order 
of magnitude, as expected from the above considerations. 

For larger volumes the agreement cannot be improved by including other contributions, 
since it is dominated by the residual truncation error of the RG-extrapolated TCSA. We 
remark that a similar evaluation of the LM series for the case of $13 operator the agreement 
of the LM series and TBA calculations gave an agreement of precision IO -13 with the TBA 
for volumes / > 7. Figure 14.21 illustrates the improvement of the LM-series by adding more 
terms, as compared to the RG-extrapolated TCSA data. 

5 A remark on finite volume form factors and /i-terms 

In this section we digress to report an issue we found when we tested the accuracy of TCSA 
by comparing numerical results for the finite volume matrix elements to predictions of the 
finite volume form factor formalism, using the methods developed in [HI E] ■ 
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Figure 4.1: Fitting the residual difference Sf 12 between the LM series evaluated up to 112 and 
the RG extrapolated TCSA result 



In [8] Pozsgay and Takacs derived a formula for finite volume form factors valid up to 
exponential corrections: 

n ,..., Jn ({Jl,...,Jn}\O(0,0)\{h,...,I m }) n _ lm . L = 



± F ° n ,...,n,i 1 i m K+™ 



H7T,01,...,tf m I 



y/pil in{.0' 1 ,..,e' n )Pi 1 ,...,i m (.0l,-,0m) 



+0 (e-» L ) 



(5.1) 



where I and J are the quantum numbers, i and j are the particle types in the given states, p 
is the density of sates and F is the infinite volume form factor, the phase factors $ have the 
form 



$ 



ii,...,i„ 



L> ■ ■ ■ j w m) 



\\ S ikh (6 k - 81 

\ k,l=l 
\ k<l 



and the ± sign corresponds to the ambiguity in choosing the branch of the square root func- 
tions. We note that the same ambiguity is present in TCSA due to the fact that the choice 
of eigenvectors is not unique. Formula (|5.ip is valid if there are no disconnected terms in the 
matrix element; for the treatment of disconnected terms cf. [9] . 

The density of states in (|5.1|) can be obtained from the finite volume quantization condi- 
tions (valid up to exponential corrections) 



Qi 1 ...i n (Oli--->0n) 
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Figure 4.2: Leclair-Mussardo series against RG-extrapolated TCSA evaluation for the VEV 
of $12 by taking account more and more contributions, in units of m\ 
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and it depends on the volume and it is positive for large volumes since its leading behavior is 
given by 



n, 



n 

fe=i 



rrii k L cosh ( 



+ 0(L 



n-l\ 



where the subleading corrections depend on the scattering phase shifts [8] . Eventually, nothing 
prevents these factors from becoming zero or negative for smaller values of the volume, which 
would therefore produce a singularity in the finite volume form factor predicted by (|5,ip . Such 
singularities have not been encountered previously; however, in the T2 model we found this 
behaviour for the state B1B2, as shown in figure l5~T1 

On the other hand, the matrix elements calculated from TCSA prove to be regular (no 
singularities jj. Therefore we should take a closer look at the finite volume form factor formula 
(|5.ip . Since it is correct only up to exponential corrections, the natural idea is that exponential 
corrections could resolve the singularities. 

In |22) Pozsgay introduced a method to describe a particular class of exponential correc- 
tions, the so-called //-terms which originate from the composite structure of the particles under 



For a non-unitary model this is not necessarily the case, as the eigenvectors can have zero norm for some 
values of the volume. These norms appear in the denominator when evaluating the matrix elements in TCSA 
[8 . This does happen in the case of the boundary Lee- Yang model when the ground state level and the first 
excited level cross 126 . 
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the bootstrap principle. Using the bound-state quantization first introduced in |27j . these can 
be described by a continuation of the quantization relations (|5.2p and the form factor formula 
(|5.ip to complex rapidities. For more details we refer the reader to the original article |22| . 
and to |28[ [29] where this method was used with great success. 

To demonstrate how the exponential corrections dissolve the singularities, let us examine 
the |{0, 0}) 12 state, which contains a B\ and a Bi particle, both with quantum number 0. The 
I?2 particle can be treated as the bound state of two B\ particles which have complex rapidity. 
We can therefore treat the state as a three- B\ state|{0, 0, 0}) m and look for a solution of (|5.2p 
with rapidities 6\ = 0, #2 = iu and #3 = —iu. The single unknown u can be obtained by 
solving 

imiLsinu + 5u(iu) + 5n(2iu) = 

and the analytic continuation of (|5.ip (including the phase factors) gives the finite volume form 
factor including the //-term contribution. Figure Iq~T1 shows the finite volume form factor with 
the naive evaluation and the correction by the //-terms against the TCSA data for the vacuum- 
B1B2 matrix elements. It can be seen clearly that the exponential corrections do resolve the 
singularities, and the agreement between the TCSA data and form factor expression is very 
good down to small volumes where other exponential corrections become relevant. 

Figure [531 shows the results for the B1B2- B\B^ diagonal matrix elemental. The diagonal 
matrix element of $1^ operator contains no singularity, since it is related to the trace of the 
energy momentum tensor, and its diagonal matrix element can be expressed using the energy 
and momentum eigenvalues of the state as computed in the approximation given by (j5.ip . 
which are finite. However even in this case including the //-terms improves the agreement 
with TCSA. For the operator $i ; 2 the conclusion is the same as in the non-diagonal case: the 
//-terms again lead to a resolution of the singularity. 

To end this section, we mention a singularity which is caused by the inclusion of //-terms, 
in a marked contrast to the above two cases where the //-terms eliminated singularities of the 
expression ()5.ip . Figure l5~3"1 shows the matrix elements multiplied by the one-particle density 
of states for the B2 particle in the vacuum- B2 case, i.e. the expression 



VM0)(0|O|{0}) 2 

where p%($$) = ITI2L 



(the extra factor y/p2(0) is put in to make the large volume asymptotics of the matrix elements 
a constant, which makes the plots easier to interpret). We see that including the //-term 
correction, while improving the agreement with TCSA, also lead to a singularity at volume 
I ~ 5.105. It is plausible that other exponential corrections would resolve these singularities. 
The next class to be taken into account are the so-called F-terms, which would lead to a change 
in the quantization condition by terms that vanish exponentially with increasing volume. For 
the energy levels, such terms can be found from the excited TBA approach of |30[ [3T| 132] ; a 
general treatment of F-terms is given in |33j . However, for the form factors the theoretical 
description of these corrections is not yet developed and is currently under active investigation 



4 Note that these matrix elements contain disconnected contributions, and therefore must be computed by 
a formula different from (|5.1[) ; for details cf. P l22ll29] . 
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Figure 5.1: Vacuum- B\B<i matrix elements for $i j2 and <J>i i2 operators in units of mi, where 
the solid line is the naive finite volume form factors result, the dashed line is the correction with 
the //-terms and the dots represent the TCSA data. All data in the plot show absolute values 
of the matrix elements to get rid of the phase ambiguities related to choice of eigenvectors in 
TCSA. 
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Figure 5.2: B1B2-B1B2 matrix elements for $i 2 and <&i j2 operators in units of m\, where 
the solid line is the naive finite volume form factors result, the dashed line is the correction 
with the /^-terms and the dots represent the TCSA data. The breaks in the lines are due to 
plotting the absolute values of the predicted matrix elements (the true value would become 
negative and the line cross below the real axis). 
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Figure 5.3: Vacuum- £>2 matrix elements multiplied by the one-particle density of states for 
the .E>2 particle for $12 and <3?i.2 operators in units of mi, where the solid line is the naive 
finite volume form factors result, the dashed line is the correction with the /i-terms and the 
dots represent the TCSA data. 
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6 Conclusions and outlook 

In this paper we investigated expectation values of local operators in integrable quantum 
field theores, in finite volume or, equivalently, at finite temperature. We considered both 
numerical and analytical approaches, and established their consistency, thereby showing that 
a combination of these methods is able to give an accurate description of these expectation 
values for all scales from small to large volumes, i.e. from the high-temperature to the the 
low-temperature regime. For detailed studies we considered the so-called T2 model, which has 
a nontrivial primary field besides the perturbing one, and all of its form factors are known 
explicitly to any number of particles. 

The ultraviolet regime can be accessed using perturbed conformal field theory, and the non- 
perturbative truncated conformal space approach extends its validity significantly and into the 
infrared regime proper. When supplemented with renormalization group improvement, it can 
be matched to methods based on the infrared date (i.e. scattering theory) with an astonishing 
precision in a wide volume/temperature range (indeed, as shown under 13.3.31 and 13.3.41 the 
RG improvement is mandatory for a proper agreement). The infrared regime can be accessed 
using the Leclair-Mussardo series, which for the perturbing operator is equivalent to the TBA 
equation. We emphasize that the matching also works for operators for which the LM series 
cannot be derived in the TBA framework, as exemplified by the $12 operator in the T2 model. 

In addition to the above main results, the paper also contains a derivation of the cut- 
off dependence of expectation values in TCSA, which has already been used before [20] . but 
without the formalism presented; the details are important in order to facilitate the use of this 
method in other calculations. Finally, while checking the accuracy of TCSA and the details of 
the T2 form factors (|B.2p using the methods developed in [8j |9] , we found interesting cases in 
which the finite volume form factor formula, which neglects corrections decaying exponentially 
with the volume, predicts singularities that are not observed in TCSA and have no theoretical 
reason to exist. It was shown how these singularities are resolved by including the class of 
exponential corrections called /i-terms. However, somewhat surprisingly, this leads to the 
appearance of new singularities, which are in turn expected to be cured by higher exponential 
corrections. These higher exponential are not yet understood and it is clear that a systematic 
description of exponential corrections to matrix elements in finite volume is necessary to settle 
all these issues. Such a description can also have other applications in understanding the 
volume/temperature dependence of physical quantities. Work is in progress in this direction. 
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A TCSA in the T 2 model 

The minimal models -M2,2n+3 have the central charge |35| 

2n(5 + 6n) 
°~ 3 + 2n 
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When perturbed by the relevant operator $13 with conformal dimensions 

l-2ra 



^1,3 = ^1,3 



3 + 2n 



for an appropriate choice of the coupling constant they give rise to the massive quantum field 

theories T n ; the n = 1 case (Ti) is the scaling Lee- Yang model. For n = 2 the central charge 

is 

68 

c = -y 

and there are two nontrivial primary field $12 and $1,3 with conformal weights hi 2 = ^1,2 = 
— 2/7 and /11.3 = fox, 3 = —3/7. We consider the theory on a Euclidean space-time cylinder of 
circumferenceL which can be mapped unto the punctured complex plane using 

2vt . _ 2tt , . 

10 = exp — (t — ix) , w = exp —(t + ix) (A.l) 

under which primary fields transform as 

2irw \ ( 2irw \ 



$(t,x)=(-^J f-^J $(w,w) 



The fields $1 j. are normalized with as 



'i,k 



(0|*i, fc (w, w)*i, fc (0, 0)|0) = — jF— (A.2) 

and the Hilbert space is given by 

n T2 = v h ® % 

fc=0,-2/7,-3/7 

where V^ (V^) denotes the irreducible representation of the left (right) Virasoro algebra with 
highest weight h. 

The Hamiltonian of the T% model takes the following form in the perturbed conformal field 
theory framework: 



H = H -X I dx$i, 3 (Q,x) 

Jo 



(A.3) 





H = 


27r ( T T c \ 

t( Lo + Lo -12) 










is the conformal Hamiltonian. When A 
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with masses mi and 
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m,2 = 2mi cos — 
5 








(A.4) 


The mass gap mi can be 


related to the 

A = 
where 


! coupling constant as 

20/7 
nm 1 

K = 0.04053795542 


HO 






(A.5) 
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and the bulk energy density is given by 

ism 



5 



Due to translational invariance of the Hamiltonian (|A.3p . the conformal Hilbert space % can 
be split into sectors characterized by the eigenvalues of the total spatial momentum 

P = ^(L -L ) 

where the operator Lq — Lq generates spatial translations and its eigenvalue is called the 
conformal spin. For a numerical evaluation of the spectrum, the Hilbert space is truncated 
by imposing a cut in the conformal energy. The truncated conformal space corresponding to 
a given truncation and fixed value s of the Lorentz spin reads 

^TCs(s,e cut ) = {|V>}G^| (Lo-Z )|V> = s|^>, ( L o + io-^)|V>} = e|V> : e < e cut } 
On this subspace, the dimensionless Hamiltonian matrix can be written as 

where energy is measured in units of the particle mass mi, I = m\L is the dimensionless 
volume parameter, 



c® = m (a.« 



1-3 

is the conformal inner product matrix and 

Bff = mi,3(w,w)\j)\ w= ^ =l (A.9) 

is the matrix element of the operator $ at the point w = w = 1 on the complex plane between 
vectors \i), \j) from %tcs( s j e cu t). The natural basis provided by the action of Virasoro 
generators is not orthonormal and therefore G^ s '~ l must be inserted to transform the left 
vectors to the dual basis. The Hilbert space and the matrix elements are constructed using 
an algorithm developed by Kausch et al. |27] , 

Diagonalizing the matrix hij we obtain the energy levels as functions of the volume, with 
energy and length measured in units of m. We considered sectors with s = 0, 1, 2, 3, and the 
maximum value of the cutoff e cu t was 30, in which case the Hilbert space contains around 
twelve thousand vectors for even, and nine thousand vectors for odd values of the conformal 
spin. Once the eigenvectors are obtained, the matrix elements of local operators can be 
computed using the exponential mapping to evaluate matrix elements; for details cf. [8]. 

B Scattering theory and form factors of the T^ model 

The form factors of T n models were constructed in pT7J using the fact that they can be ob- 
tained as reductions of sine-Gordon theory at a particular value of the coupling [38J. In these 
particular models, only the breather states are retained in the spectrum. For T<i only the first 
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two breathers B\ and B<i are retained in the spectrum, and their two-particle S matrices take 
the form 



Sn(0) 
with 



Si 2 (0) 



{*}« 



11 (S 
5J J5 

sinh 9 + % sin irx 
sinh 9 — i sin -zra; 



S22OT 



-V - 

5 1.15 
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The form factors of primary fields with the fundamental particle can be obtained from the 
form factors of exponential operators in sine-Gordon theory with the first breather. The latter 
are given by 
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where the parameter £ is 
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giving the minimal B\B\ form factoiO, while cr\' denotes the elementary symmetric polyno- 
mial of n variables and order k defined by 

n n 



t=l 



fe=0 



5 The formula for the function v is in fact independent of JV; choosing JV large extends the width of the 
strip where the integral converges and also speeds up convergence. 
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and Q a {P) is the vacuum expectation value of the field which is known exactly [39]. Form 
factors of higher breathers can be computed by representing them as bound states of B\ 
particles; a useful formula for their evaluation can be found in Appendix A of [40] . 
The models T n correspond to restriction at the coupling 



2n + l 



For the T<i model, the restricted spectrum is composed of the first and the second breathers 
B\ and Bi- The form factors of the operators $12 and $^3 can be obtained as the cases 
a = 1 and a = 2 from formula (|B.2p [37J; however, the vacuum expectation value needs to be 
replaced by the exact vacuum expectation value of the minimal model fields, derived in |41) : 



($1,2) 
($1,3) 



-2.3251365527- 



x im 



-4/7 
1 



2.2695506880 • • • X m 



-6/7 
1 



(B.4) 



A form factor containing n B\ particles and 771 B2 particles can be evaluated using a funda- 
mental form factor (JB.2J) containing n + 2m B\ particles; therefore the connected form factor 
in (|2.3p can be evaluated from a 2n + 4m-particle fundamental form factor. 



References 

A. Leclair and G. Mussardo, "Finite temperature correlation functions in integrable 



[1 
[2 

[3] 
[4] 
[5] 
[6] 
[7 



[9 



QFT," Nucl.Phys. B552 (1999) 624-642 arXiv:hep-th/9902075 [hep-th] 



C.-N. Yang and C. Yang, "Thermodynamics of one-dimensional system of bosons with 
repulsive delta function interaction," J.Math.Phys. 10 (1969) 1115-1122. 

A. Zamolodchikov, "Thermodynamic Bethe Ansatz in Relativistic Models. Scaling 



Three State Potts and Lee- Yang Models," Nucl.Phys. B342 (1990) 695-720 



M. Karowski and P. Weisz, "Exact Form-Factors in (l+l)-Dimensional Field Theoretic 
Models with Soliton Behavior," \Nucl. Phys. B139 (1978J455J 

A. N. Kirillov and F. A. Smirnov, "A representation of the current algebra connected 



with the SU(2) invariant Thirring model," Phys. Lett. B198 (1987) 506-510 



F. A. Smirnov, "Form-factors in completely integrable models of quantum field theory," 
Adv. Ser. Math. Phys. 14 (1992) 1-208. 

H. Saleur, "A comment on finite temperature correlations in integrable QFT," 



Nucl. Phys. B567 (2000) 602-610, arXiv:hep-th/9909019 



B. Pozsgay and G. Takacs, "Form factors in finite volume I: form factor bootstrap and 



truncated conformal space," Nucl. Phys. B788 (2008) 167-208, 



arXiv : 0706 . 1445 [hep-th] 



B. Pozsgay and G. Takacs, "Form factors in finite volume. II. Disconnected terms and 



finite temperature correlators," Nucl.Phys. B788 (2008) 209-251, 



|arXiv: 0706 .3605 [hep-th] 



20 



[10] 

[11] 

[12] 

[13] 
[14] 

[15] 
[16] 
[17] 

[18] 
[19] 

[20] 

[21] 
[22] 

[23] 
[24] 



B. Pozsgay, "Mean values of local operators in highly excited Bethe states," 



J.Stat.Mech. 1101 (2011) P01011, arXiv: 1009.4662 [hep-th] 



M. Kormos, A. Trombettoni, and G. Mussardo. "Expectation Values in the Lieb-Liniger 
Bose Gas," \Phys. Rev. Lett. 103 (2 009) 210404. Selected for the December 2009 issue of 
Virtual Journal of Atomic Quantum Fluids, Vol.1, Issue 6. 

M. Kormos, G. Mussardo, and A. Trombettoni, "Local Correlations in the Super 
Tonks-Girardeau Gas," Phys.Rev. A83 (2011) 013617, 
arXiv : 1008 . 4383 [cond-mat . quant-gas] , 

B. Pozsgay, "Local correlations in the ID Bose gas from a scaling limit of the XXZ 
chain," \J.Stat.Mech. 1111 (2011) P11017[ |arXiv : 1108 . 6224 [cond-mat . stat-mech] [ 

F. H. L. Essler and R. M. Konik, "Finite-temperature dynamical correlations in massive 
integrable quantum field theories," [7. Stat. Mech. 090 9 (2009) P09018, 
larXiv : 0907 . 0779 [cond-mat . str-el] [ 

F. H. L. Essler and R. M. Konik, "Finite-temperature lineshapes in gapped quantum 



spin chains," \Phys. Rev? B78 (2008) 100403, [arXiv : 071 1 ■ 2524 [cond-mat . str-el]] 
B. Pozsgay and G. Takacs, "Form factor expansion for thermal correlators," 



J.Stat.Mech. 1011 (2010) P11012, arXiv: 1008.3810 [hep-th] 



I. M. Szecsenyi and G. Takacs, "Spectral expansion for finite temperature two-point 



functions and clustering," J.Stat.Mech. 1212 (2012) P12002 



arXiv: 1210.0331 [hep-th] 



V. P. Yurov and A. B. Zamolodchikov, "Truncated conformal space approach to scaling 



Lee- Yang model," Int. J. Mod. Phys. A5 (1990) 3221-3246. 



G. Feverati, K. Graham, P. A. Pearce, G. Z. Toth, and G. M. T. Watts, "A 
renormalization group for the truncated conformal space approach," 



J. Stat. Mech. 2008 no. 03, (2008) P03011, arXiv:hep-th/0612203 [hep-th] 



R. M. Konik and Y. Adamov, "A Numerical Renormalization Group for Continuum 
One-Dimensional Systems," \Phys. Rev. Lett. 98 (2007) 147205[ 
arXiv : cond-mat/0701605 [cond-mat . str-el] . 

P. Giokas and G. Watts, "The renormalisation group for the truncated conformal space 
approach on the cylinder," arXiv: 1106.2448 [hep-th] , 

B. Pozsgay, "Luscher's mu-term and finite volume bootstrap principle for scattering 
states and form factors," \Nucl.Phys. B802 (2 008) 43 5-457[ 



arXiv : 0803 . 4445 [hep-th] 



T. Palmai and G. Takacs, "Diagonal multi-soliton matrix elements in finite volume," 
Phys. Rev. D 87 (2013) 045010, |arXiv: 1209.6034 [hep-th]] 



A. Zamolodchikov, "Integrals of Motion and S Matrix of the (Scaled) T=T(c) Ising 



Model with Magnetic Field," Int.J.Mod.Phys. A4 (1989) 4235 



27 



[25] 

[26] 

[27] 

[28] 
[29] 
[30] 

[31] 
[32] 
[33] 

[34] 
[35] 

[36] 

[37] 
[38] 
[39] 

[40] 



G. Feverati, F. Ravanini, and G. Takacs, "Truncated conformal space at c = 1, 
nonlinear integral equation and quantization rules for multi - soliton states," 
\Phys.Let t. B430 (1998) 264-273[|arXiv:hep-th/ 9803104 [hep-th]] 

P. Dorey, M. Pillin, R. Tateo, and G. Watts, "One point functions in perturbed 
boundary conformal field theories," Nucl.Phys. B594 (2001) 625-659, 
|arXiv : hep-th/0 007077 [hep-th] | 

H. Kausch, G. Takacs, and G. Watts, "On the relation between Phi(l,2) and Phi(l,5) 



perturbed minimal models," Nucl.Phys. B489 (1997) 557-579, 



arXiv:hep-th/9605104 [hep-th] 



G. Feher and G. Takacs, "Sine-Gordon form factors in finite volume," 



Nucl.Phys. B852 (2011) 441-467, arXiv: 1106. 1901 [hep-th] 



G. Takacs, "Determining matrix elements and resonance widths from finite volume: the 



dangerous //-terms," JE EP 1111 (2011) 113] |arXiv: 1110.2181 [hep-th] 

V. V. Bazhanov, S. L. Lukyanov, and A. B. Zamolodchikov, "Integrable quantum field 



theories in finite volume: Excited state energies," Nucl.Phys. B489 (1997) 487-531 
|arXiv : hep^ th/9607099 [hep-th] | 

P. Dorey and R. Tateo, "Excited states by analytic continuation of TBA equations," 
\Nucl.Phy s. B482 (1996) 639-659[rarXiv:hep-th7 960716 7 [hep-th] | 

P. Dorey and R. Tateo, "Excited states in some simple perturbed conformal field 
theories," \Nucl.Phys. B515 (1998) 575-623[ |arXiv:hep-th/9706140 [hep-th]] 

Z. Bajnok and R. A. Janik, "Four-loop perturbative Konishi from strings and finite size 
effects for multiparticle states," \Nucl.Phys. B807 (2009) 625-650[ 
|arXiv: 0807 .0399 [hep-th] | 

B. Pozsgay, in preparation. 

A. Belavin, A. M. Polyakov, and A. Zamolodchikov, "Infinite Conformal Symmetry in 
Two-Dimensional Quantum Field Theory," \Nucl.Phys. B241 (1984) 333-380[ 

V. Fateev, "The Exact relations between the coupling constants and the masses of 
particles for the integrable perturbed conformal field theories," 



Phys.Lett. B324 (1994) 45-51 



A. Koubek, "Form-factor bootstrap and the operator content of perturbed minimal 



models," Nucl.Phys. B428 (1994) 655-680 arXiv:hep-th/9405014 [hep-th] 



N. Reshetikhin and F. Smirnov, "Hidden Quantum Group Symmetry and Integrable 
Perturbations of Conformal Field Theories," Commun.Math.Phys. 131 (1990) 157-178 



S. L. Lukyanov and A. B. Zamolodchikov, "Exact expectation values of local fields in 



quantum sine-Gordon model," Nucl.Phys. B493 (1997) 571-587 



|arXiv:hep-th/9611238 [hep-th] 



G. Takacs, "Form factor perturbation theory from finite volume," 



Nucl. Phys. B825 (2010) 466-481[|arXiv: 0907 .2109 [hep-th]] 



28 



[41] V. Fateev, S. L. Lukyanov, A. B. Zamolodchikov, and A. B. Zamolodchikov, 

"Expectation values of local fields in Bullough-Dodd model and integrable perturbed 
conformal field theories," \Nucl.Phys. B516 (1998) 652-674[ 



|arXiv : hep"^t h/9709034 [he p-th] 



29 



